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Abstract 

Let s = (si, . . . , s m ) and t = (ti, . . . ,t n ) be vectors of nonnegative integer- valued 
functions of m, n with equal sum S = Y^=i s « = Sj=i Let ^( s > *) be the number 
ofmxn matrices with nonnegative integer entries such that the ith row has row 
sum S{ and the jth column has column sum tj for all i,j. Such matrices occur in 
many different settings, an important example being the contingency tables (also 
called frequency tables) important in statistics. Define s = maxj Sj and t = maxj tj. 
Previous work has established the asymptotic value of M(s, i) as m, n — > oo with s 
and t bounded (various authors independently, 1971-1974), and when all entries of 
s equal s, all entries of t equal t, and m/n,n/m, s/n > cj log n for sufficiently large c 
(Canfield and McKay, 2007). In this paper we extend the sparse range to the case 
st = o(S 2 / 3 ). The proof in part follows a previous asymptotic enumeration of 0-1 
matrices under the same conditions (Greenhill, McKay and Wang, 2006). We also 
generalise the enumeration to matrices over any subset of the nonnegative integers 
that includes and 1. 



1 Introduction 

Let s = s(m, n) = (si, . . . , s m ) and t = t(m, n) = (t±, . . . , t n ) be vectors of nonnegative 
integers with equal sum S = YlT=i s i = Y^=i^j- Let M(s, t) be the set of all m x n 
matrices with nonnegative integer entries such that the ith row has row sum and the 
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jth column has column sum tj for each Then define M(s,t) = \Ai(s,t)\ to be the 
number of such matrices. 

Our task in this paper is to determine the asymptotic value of M(s,t) as m, n — > oo 
under suitable conditions on s and t. 

The matrices Ai(s,t) appear in many combinatorial contexts; see Stanley [T5J Chap- 
ter 1] for a brief history. A large body of statistical literature is devoted to them under 
the name of contingency tables or frequency tables; see [HI E| for a partial survey. In theo- 
retical computer science there has been interest in efficient algorithms for the problem of 
generating contingency tables with prescribed margins at random, and for approximately 
counting these tables. See for example [TJ [8j [T3j . 

The history of the enumeration problem for nonnegative integer matrices is surveyed 
in [5], while a history for the corresponding problem for 0-1 matrices is given in [TT]. Here 
we recall only the few previous exact results on asymptotic enumeration for nonnegative 
integer matrices. Define s = maxj Si and t = maxj tj. The first non-trivial case s± = ■ ■ ■ = 
s m = h = ■ ■ ■ = t n = 3 was solved by Read [14] in 1958. During the period 1971-74, 
this was generalised to bounded s, t by three independent groups: Bekessy, Bekessy and 
Komlos [2], Bender [3], and Everett and Stein [9], under slightly different conditions. 

In the case of denser matrices, the only precise asymptotics were found by Canfield 
and McKay [5] in the case that the row sums are all the same and the column sums are all 
the same. Let M(m, s; n, t) = M((s, s, . . . , s), (t, t, . . . , t)), where the vectors have length 
m, n, respectively, and ms = nt. 

Theorem 1.1 (^J Theorem 1]). Let s = s(m, n), t = t(m, n) be positive integers satisfying 
ms = nt. Define X = s/n = t/m. Let a,b > be constants such that a + b < |. Suppose 
that m, n — > oo in such a way that 



Then A(m, s; n, t) = 0(n b )(m + n) as m, n — > oo. □ 

Canfield and McKay conjectured that in fact < A(m, s;n,t) < 2 for all s,t > 1. 
The results in the present paper establish that conjecture for sufficiently large m, n in the 
case st = o{{mn) l ' b ). (See Corollary Q) 




(1.1) 



Define Aim, s; n, t) by 
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The main result in this paper is the asymptotic value of M(s, t) for st = o(S 2 ^ 3 ). Our 
proof uses the method of switchings in a number of different ways. In several aspects our 
approach is parallel to that which provided our previous asymptotic estimate of N(s,t), 
the number of 0-1 matrices in the class A4(s, t). We now restate that result for conve- 
nience. For any x, define [x]o = 1 and for a positive integer k, [x]k = x(x—l) ■ ■ ■ (x—k+1). 
Also define 

m n 

Sk = 2^ l S i\ki Tk = [tj]k 

t=l j=l 

for k > 1. Note that Si = T x = S. 

Theorem 1.2 ([TTJ Corollary 5.1]). Let s = s(m,n) = (si, ...,s m ) and t = t(m,n) = 
(ti,...,t n ) be vectors of nonnegative integers with equal sum S = Y^ILi s « = Sj=i*j- 
Suppose that m, n — > oo, S — > oo and 1 < st = o(S 2 ^ 3 ). Then 

Ar , ,x S\ I S2T2 S2T2 S3T3 S2T2{S2 + T2) 

N ( s ^) = TWi , T-rn 71 ex P ?To9~ - TTm" + 



n'l^in^i^'! V 2S 2 2S 3 3S 3 4S 4 

S 2 T 3 + S 3 T 2 , S 2 T 2 . ^ S 3 t 3 



2^ + is? + °[-* ]l D 



We now state our main result, which is the asymptotic value of M(s, t) for sufficiently 
sparse matrices. Note that the answer is obtained by multiplying the expression for 
N(s,t) from Theorem 11.21 by a simple adjustment factor. 

Theorem 1.3. Let s = s(m, n) = (si, . . . , s m ) and t = t(m, n) = (ti, . . . , t n ) be vectors 
of nonnegative integers with equal sum S = Y^Li s i = Sj=i tj- Suppose that m,n — ► oo ; 
S -> 00 and 1 < st = o(S 2 / 3 ). Then 

S0T0 S,U „ /s 3 t 3 



M(s, t) = iV( S , t) exp + ^ + °( 



S 2 



S\ ( S2T2 S2T2 S3T3 S^I^nS^ + T2 



n™iSi!nj=i^' ! \2S 2 2S 3 3S 3 4S 4 



S 2 2 r 3 + S 3 T| | S 2 2 T| | g ^ 3 t 3 



2S 4 2S 5 V S 

Proof. The proof of this theorem is presented in Sections [2] and [31 First we show that 
the set of matrices in A4(s, t) with an entry greater than 3 forms a vanishingly small 
proportion of M.(s,i). We also show that it is very unusual for an element of A4(s,t) 
to have a "large" number of entries equal to 2 or a "large" number of entries equal to 3, 
where "largeness" is defined in Section [2j We establish these facts using switchings on 
the matrix entries. This allows us to concentrate on matrices in A4(s,t) with no entries 
greater than 3 and not very many entries equal to 2 or 3. 
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We then proceed in Section [3]to compare the number of these matrices with the number 
N(s, t) of {0, l}-matrices with row sums s and column sums t. We do this by adapting 
the results from [TT] used to prove Theorem 11.21 These calculations are carried out in 
the pairing model, which is described in Section [3j Our theorem follows on combining 
Lemmas 13.11 13.21 and Corollary 13.81 □ 

In the semiregular case where s, = s for 1 < i < m and tj = t for 1 < j ' < n, 
Theorem 11.31 says the following. 

Corollary 1.4. Suppose that m, n — > oo and that sm = tn = S for nonnegative integer 
functions s = s(m,n), t = t(m, n) and S = S(m,n). If 1 < st = o(S 2 ^ 3 ) then 

M{m, s; n, t) = 

S\ / ( a -i)(t-i) (s-l)(t-l)(2st-s-t-10) | c (s 3 t 3 \\ n 



(s\) m (t\) n L \ 2 12S \S 2 

For some applications the statement of Theorem 11.31 is not very convenient. In Sec- 
tion |4] we will derive an alternative formulation, very similar to one given for N(s, t) 
in [TT]. For k — 2,3, define 



= nr..... , n\ > A S i ~ 5 / m ) 



S(mn + S) 



i=i 

n 

ran 



S(mn + S) 

3=1 



To motivate the definitions, recall that S/m is the mean value of Sj and S/n is the mean 
value of tj, so these are scaled central moments. We will prove Corollary 14.11 stated in 
Section HI which has the following special case. 

Corollary 1.5. Under the conditions of Theorem ] 1.!% if (1 + /^(l + v 2 ) = 0(5 1 / 3 ) then 

n (" + r n (T 1 ) , 

M(s,t) = -= exp(i(l - *)(! - ft,) +o(^)). □ 

Corollary 11.51 has an instructive interpretation. Following [5], we write M(s,t) = 
MP X P 2 E, where 

mn+ s s ~ 1 ), p^m-^^j- 1 ), ft=^ i n( mr ^ i " 1 )' 

£ = exp - A 2 )(l - ^2) + O ) . 
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Clearly, M is the number ofmxn nonnegative matrices whose entries sum to S. In the 
uniform probability space on these M matrices, Pi is the probability of the event that the 
row sums are given by s and P 2 is the probability of the event that the column sums are 
given by t. The final quantity E is thus a correction to account for the non-independence 
of these two events. 

Finally, in Section Owe show how to generalise Theorem II .31 to matrices whose entries 
are restricted to any subset of the natural numbers that includes and 1. 

A note on our usage of the 0( ) notation in the following is in order. Given a fixed 
function f(S) = o(S 2 ^ 3 ), and any quantity <p that depends on any of our variables, 0(0) 
denotes any quantity whose absolute value is bounded above by \aj>\ for some constant c 
that depends on / and nothing else, provided that 1 < st < f(S). 



2 Switchings on matrices 

In this section we will show that the condition st = o(S 2 ^ 3 ) implies that most matrices 
have no entries greater than 3. We also find bounds on the number of entries equal to 2 
or 3. Our tool will be the method of switchings, which we will analyse using results of 
Fack and McKay [10] from which we will distill the following special case. 

Theorem 2.1. Let G = (V, E) be a finite simple acyclic directed graph, with each v G V 
being associated with a finite set C(v), these sets being disjoint. Suppose that S is a 
multiset of ordered pairs such that for each (Q, R) G S there is an edge vw G E with 
Q G C(y) and R G C(w). Further suppose that a, b : V — > R are positive functions such 
that, for each v G V , 

\{(Q,R)eS\QeC(v)}\>a(v)\C(v)\, 

2.1 

\{(Q,R)eS\ReC(v)}\ <b(v)\C(v)\, 

where the left hand sides are multiset cardinalities. Let 7^ Y C V. Then there is a 
directed path Vi, v%, . . . , in G, where V\ G Y and t> fc is a sink, such that 

E v ev\C(v)\ ~ Ei<><kNW ' 

where N(vi) is defined by 

jv(«,) = 1, 

%i = V\"ti p < < < *). 

b[v 2 ) ■ ■ ■ b{vi) 

Proof. This follows from Theorems 1 and 2 of [TD]. □ 



5 



For D > 2, a D -switching is described by the sequence 

(Q; (io, Jo), (iiji),---, fa, 3d)) 

where Q is a matrix in A4(s,t) and (i , j ), (ii,ji), . . . , (iD,jo) is a (D+l)-tuple of posi- 
tions such that 

• the rows i , . . . ,i D are all distinct and the columns j , ■ ■ ■ ,3d are all distinct; 

• there is a D in position (i , j ) of Q; 

• the entries in positions (ie,je) of Q are not equal to or D + 1, for 1 < £ < D; 

• there is a in position (ie,jo) and position (io,je) of Q for 1 < £ < D. 

This D-switching transforms Q into a matrix _R G (s, t) by acting on the (D+l) x (-D+1) 
submatrix consisting of rows (io, . . . ,io) and columns (jo, • • • , 3d) as follows: 



/ 



Q = 



D 
qi 
q 2 







\ 



V 



/ 



1 1 

1 qi-1 

1 92-1 



\ 



V 



qo- 



= R. 



I 



Matrix entries not shown can have any values and are unchanged by the switching oper- 
ation. Notice that the D-switching preserves all row and column sums and reduces the 
number of entries equal to D by at least 1 and at most D + l. The number of entries 
greater than D is unchanged. 

A reverse D -switching, which undoes a D-switching (and vice- versa), is described by 
a sequence (R; (io, Jo), • • • , («d, jo)) where R e M(s,t) and (i ,jo), (hji), • • • , (id, 3d) 
is a (D+l)-tuple of positions such that 

• the rows i , . . . ,i D are all distinct and the columns jo, ■ ■ ■ ,3d are all distinct; 

• there is a zero in position (i , j ) of R; 

• the entries in positions (ie,je) of R are not equal to D, for 1 < £ < D; 

• there is a 1 in position (ig,jo) and position (io,je) of R for 1 < £ < D. 
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Lemma 2.2. Let D > 2 and let Q e Ai(s, t) have at least K > 2st non-zero entries 
that are not greater than D, and at least J entries equal to D. Then there are at least 
J(K — 2st) D D-switchings and at most SdTjj reverse D-switchings that apply to Q. 

Proof. First consider D-switchings. We want a lower bound on the number of (D+l)- 
tuples (io, Jo) j • • • j (1d> Jd) of indices where a D-switching may be performed. There are at 
least J ways to choose the position (20, jo)- Then we can choose the remaining positions 
one at a time, avoiding choices which violate the rules. The choice of the last position 
(^d, Jd) is the most restricted, so we bound that. By assumption, there are at least K 
nonzero entries in Q that are not greater than D. Of these we must exclude the entry in 
position (io, jo) as we ll as entries in the same column as a nonzero entry in row %q other 
than column jo (at most (s — D)t positions), entries in the same row as a nonzero entry in 
column jo other than row io (at most (t — D)s positions), and entries in row ii or column 
je for 1 < i < D — 1 (at most (D — l)(s+t — 2) positions). Overall, we can choose position 
(in, Jd) in at least 

K - 1 - (s - D)t - (t - D)s - (D - l)(s + t - 2) > K - 2st 

ways, and as we noted this also applies to each of the less restricted positions (ie,je), 
where 1 < i < D. Hence at most JiK — 2st) D D-switchings involve Q. 

Next consider reverse D-switchings. An ordered sequence of D entries in the same row 
which equal 1 may be chosen in at most Su ways, and an ordered sequence of D entries 
in the same column which equal 1 may be chosen in at most Tjj ways. Some of these 
choices will not give a legal position for a reverse D-switching, but SdT d is certainly an 
upper bound. □ 

Our first application of switchings will be to show that only a vanishing fraction of 
our matrices have any entries greater than 3. For j > and D > 2, let M.d{]) be the 
set of all matrices in A4(s, t) with exactly j entries equal to D and none greater than D. 
Define M D (>0) = \Jj >0 M D (j), and note that M D+1 (0) = M D (0) U M D (>0). 

Lemma 2.3. Suppose that 1 < st — o(S 2 ^ 3 ). Let U\ = Ui(s,t) be the set of all matrices 
in Ai(s,t) which contain an entry greater than 3. Then 

\U l \/M(s,t) = 0{s i t i /S 2 ). 

Proof. The largest possible entry of a matrix in Ai(s,t) is A = min{s,t}. We will apply 
Theorem 12.11 to successively bound the possibility that the maximum entry is D, for 
D = A,A-1,...,4. 

Fix D with 4 < D < A. Define a directed graph G = (V, E) with vertex set V = 
{v , v i, v 2, . . . } and edge set E = {vjVi \ j — D — 1 < i < j — 1}. Associate each Vi with the 
set C{v i) = M. z?(i). Define S to be the set of pairs (Q, R) related by a D-switching, where 
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Q G Vj,R G Vi for some VjVi G E. Define Y = {vi,v 2 , ■ ■ ■} C V. Note that SdT d > 
since -D < A. 

We can now use Theorem 12.11 to bound 

\M D {>0)\ _ Ev&\C(v)\ 



\M D+1 (0)\ EvevlCivW 

once we have found positive functions a, b : V — > K satisfying (12.11) . These are provided 
by Lemma 12.21 with J = j and K = S/D, the latter being clear since there are no 
entries greater than D and the total of all the entries is S. We have S/D > 2st since 
D < A < (st) 1 / 2 . Thus we can take a( Vj ) = j{S/D - 2st) D and b( Vj ) = S D T D . 

Theorem 12. II tells us that, unless A^d(>0) = 0, there is a directed path v tl , v t2 , . . . , v tq , 
where q > 1 and t\ > t 2 > ■ ■ ■ > t q — (since v is the only sink) such that (12.21) holds. 

Hence, using the values of iV as given in Theorem 12.11 we have 

\M D (>0)\ < N(v tq _ 1 ) + --- + N(v tl ) 



\M D+1 (0)\ ^ N(v tq ) + --- + N(v t2 ) 
< max ——, — — 

l<i<g N{V ti ) 

b{M D {U)) 



max 



< 



i<i<? a(M D (ti-i)) 

SdTd 
(S/D-2st) D ' 



Let £d denote this upper bound: that is, £d = SdTd/(S/D — 2st) D for 4 < D < A. Note 
that ^ 4 = 0(sH 3 /S 2 ). For 4 < D < A, we have £ D > and 

eo+i <ct f(D + l) D+1 \ (S-2stD) D 



,D 



D D ) {S - 2st{D + 1)) D +! 



S - 2st(D + 1) V 5 - 2stD J 
= o(l) 

uniformly over D, where the last step uses the observation that A < (st) 1 / 2 = o(S 1 ^ 3 ). 

Since U x = A^ 4 (>0)UA^ 5 (>0)U- • -UA^^(>0) andA^ D+ i(0) C M(s,t) for 4 < D < A, 
we have \U 1 \/M(s, t) < + ^ 5 + • • • + £ A = 0(sH 3 /S 2 ) as required. □ 

We may therefore restrict our attention to matrices with no entry greater than 3. 
Next we find upper bounds on the numbers of entries equal to 2 or 3 which hold with 
high probability. 
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Define 



22 if S 2 T 2 < S 7 ' A , 

N 2 = { [log S] if S 7 ^ < S 2 T 2 <skoS 2 log S, 

\5600S 2 T 2 /S 2 ] if S 2 log S < S 2 T 2 ; 
N 3 = max([logS], \ 230000 S 3 T 3 / S 3 ]) . 

(Here and throughout the paper we have not attempted to optimise constants.) 
We will use the following lemma. 

Lemma 2.4. Let k be a positive integer and let q and n be positive real numbers such 
that n > kq. Then 

n{n — q) ■ ■ ■ {n — (k — l)q) > (n/e) k . 

Proof. Dividing the left side by n k gives, for n > kq, 

fc-i / fc-i 



K — -L / K—X \ 

-iq/n) = expf ^log(l - iq/n) J 

(jf 



i=0 i=0 

"k 

>expl / log(l — xq/n) dx 



= exp(— k — (n/q — k) log(l — kq/n)) 
> exp(— k). 

The second line holds because log(l — xq/n) is a decreasing function for x G [0, k]. The 
case n = kq follows by continuity. □ 

Lemma 2.5. Let 1 < st = o(S 2 ^ 3 ). Then, with probability 1 — 0(s 3 t 3 /S 2 ), a random 
element of Ai(s,t) has no entry greater than 3, at most N% entries equal to 3, and at 
most N 2 entries equal to 2. 

Proof. In view of Lemma [273| we may restrict our attention to the set Ai^O) of all matrices 
in Ai(s,t) with maximum entry at most 3. We will start by applying 3-switchings as in 
Lemma [2.31 but the analysis will be more delicate. 

In applying Theorem 12.11 we have V = {v ,Vi, . . . }, with Vh associated with A4 3 (h), 
and Y = {vh \ h > N 3 }. For sufficiently large S, we have from Lemma [2.21 that we can 
take a(vh) = ^§hS 3 and b(vh) = S^T^. If S3T3 = then entries equal to 3 are impossible, 
so we assume that S 3 T 3 > 0. Define <p = 28S 3 T 3 / S 3 . 

According to Theorem 12 .![ there is a sequence 

hi > h 2 > ■ ■ • > h q = 0, 
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with hi > N 3 and h^% — 4 < hi < hi-i for all i, such that 

\M 3 (>N 3 )\ \M 3 (>N 3 )\ 
\M{a,t)\ ~ \M,(0)\ 

N(h e ) + Nih-x) + ■ ■ ■ + Nih] 



< 



N(h q ) + N(h g _ 1 ) + --- + N(h 



where £ is the largest index such that h^ > N 3 + 1 and N(hi) = hi ■ ■ ■ hi-\ip 1+1 for all i. 
Define u = |_z 1°§ $\ ■ Since N 3 > [log S] , we have £ + u < q. Also, for < i < £ — 1, 



Therefore, 



N(h e + u -i) ~ N(h i+U ) 



|A4 3 (>iV 3 )| < N(ht) + iV(A^-i) + • • • + iV(/ii) 



" N{h i+U ) + N{h i+U - 1 ) + ---N{h u+1 ) 
N{h e ) 



< 



N(h e+U ) 



< 



heht + i ■ ■ ■ he +u -\ 



(N 3 + l)(N 3 -3)---(N 3 -Au + 5Y 
Since iV 3 + 1 > 4w we can apply Lemma 12.41 to obtain the bound 

\M 3 (>N 3 )\ < ( <pe 



\M(s,t)\ ~ \N 3 + 1, 
Now iV 3 > 230000S 3 T 3 /S 3 > 8214 and u > \ \ogS - 1, so this upper bound is at most 



8214 



7 log 5-1 



This shows that with probability 0(s 3 t 3 / S 2 ) there are at most N 3 entries equal to 3, as 
required. 

To bound the number of entries equal to 2, we proceed in the same manner using 
2-switchings, working under the assumption that there are at most iV 3 entries equal to 3 
and none greater than 3. In applying Lemma [2.2[ we can take K — \{S — 3iV 3 ), so that 
(K — 2st) 2 > |5* 2 for sufficiently large S. Define ip = 5S 2 T 2 /S 2 . Arguing as above we 
find a sequence 

di > d 2 > ■ ■ ■ > d r = 0, 

with the following properties: (i) d\ > N 2 and d{-i — 3 < di < di-i for all i, and (ii) if p is 
the greatest integer such that d p > N 2 then, for any w with < w < r — p, the probability 
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that there are more than N 2 entries equal to 2, subject to there being at most N 3 equal 
to 3, is bounded above by 

(2.3) 



First suppose that S 2 T 2 < S 7/i , so that N 2 = 22 and ^ < 5S~ 1/4 . Since d p > N 2 + 1 = 
23, it follows that r — p > 8. Taking w = 8 in (12.31) gives 

\ Y \ ^ -C(S-*). 

\M(s,t)\ d p d p+1 ■ ■ ■ d p+7 

Now suppose that S 2 T 2 > 5* 7//4 . Then N 2 > [log 5*] so we can take w = log .SJ . 
Arguing as above by applying Lemma [2^1 to (12.31) . we obtain the bound 0(S~ 2 ) again. 
This completes the proof. □ 

From now on we proceed in two cases, as in [TT]. Say that the pair (S 2 , T 2 ) is substantial 
if the following conditions hold: 

• 1 < st = o(S 2 / 3 ), 

• S 2 > s log 2 S and T 2 > t log 2 S, 

• S 2 T 2 > (stf/ 2 S. 

Lemma 2.6. If 1 < st = o(S 2 ^ 3 ) and (S 2 , T 2 ) is insubstantial, then with probability 
1 — 0(s 3 t 3 / S 2 ) , a random element of Ai(s,t) has no entry greater than 3, at most one 
entry equal to 3 and at most two entries equal to 2. 

Proof. The absence of entries greater than 3 follows from Lemma 12. 3[ We can also, by 
Lemma T2.51 assume that the number of entries equal to 2 or 3 is o(S). Therefore, most of 
the matrix entries are or 1. Let M be the set of all matrices in Ai(s, t) with no entries 
greater than 3, at most N 2 entries equal to 2 and at most N 3 entries equal to 3. 

To bound the number of entries equal to 2 or 3 even more tightly, as this lemma 
requires, we employ D-switchings (D = 2, 3) with the additional restriction that q\ = 
■ • • = Id — 1- This ensures that these restricted D-switchings reduce the number of 
entries equal to D by exactly one and do not create any new entries equal to 2 or 3. 

Let N"(h) be the number of matrices in Af with h entries equal to 3. If Q is such a 
matrix then the number of restricted 3-switchings applicable to Q is hS 3 (l + o(l)) and the 
number of reverse restricted 3-switchings is at most S3 T 3 . (This follows using arguments 
similar to those in Lemma [2. 2 \ since there are S — o(S) entries equal to 1.) Therefore, if 
the denominator is nonzero, 

N " {h) -O^m. (2.4) 



N"(h-1) v ' hS 3 
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We can now easily check that each of the three causes of insubstantiality (namely, S 2 < 
slog 2 S, T 2 < thg 2 S, and S 2 T 2 < (st) 3 / 2 S) imply that 



0(s 3/2 t 3/2 /S) = o(l). 



S 3 

Hence ( 12.41) implies that 

^ h - 0(s 3 t 3 /S 2 ) 



N"(0) 

In precisely the same way, using restricted 2-switchings, we find that 

N'(0) ~ { 7 h 

where N'(d) is the number of matrices in Af with d entries equal to 2 and at most one 
entry equal to 3. The lemma follows. □ 



3 From pairings to matrices 

The remainder of the paper will involve calculations in the pairing model, which we now 
describe. (This model is standard for working with random bipartite graphs of fixed 
degrees: see for example [12].) Consider a set of S points arranged in cells X\, x 2 , ■ ■ ■ , x m , 
where cell x^ has size Sj for 1 < % < m, and another set of S points arranged in cells 
yi, y 2 , . . . , y n where cell yj has size tj for 1 < j < n. Take a partition P (called a pairing) 
of the 2S points into S pairs with each pair having the form (x, y) where x G Xi and 
y G yj for some i,j. The set of all such pairings, of which there are S\, will be denoted 
by V(s,t). We work in the uniform probability space on V(s,t). 

Two pairs are called parallel if they involve the same two cells. A parallel class is a 
maximal set of mutually parallel pairs. The multiplicity of a parallel class (and of the 
pairs in the class) is the cardinality of the class. As important special simple pair 

is a parallel class of multiplicity one, a double pair is a parallel class of multiplicity two, 
while a triple pair is a parallel class of multiplicity three. 

Each pairing P G V(s, t) gives rise to a matrix in Ai(s, t) by letting the (i, j)-th entry 
of the matrix equal the multiplicity of the parallel class from x« to yj in P. 

In [UJ we noted that the number of pairings which gives rise to each 0-1 matrix in 
A4(s,t) depends only on s and t and is independent of the structure of the matrix. Hence 
the task of counting such matrices reduces to finding the fraction of pairings that have no 
multiplicities greater than 1. 

More generally, matrices in Ai(s, t) correspond to different numbers of pairings. For a 
pairing P G V(s, t), define the multiplicity vector of P to be a(P) = {a 2 , a 3 , . . . ) where a r 
is the number of parallel classes of multiplicity r. Also define the weight of P as 

w(P) = (2!) aa (3!) a3 (4!) a4 • ■ • 
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For Q G A4(s,t), define w(Q) and a{Q) to be the common weight and multiplicity vectors 
of the pairings that yield Q. 

By elementary counting, a matrix Q £ A4(s, t) corresponds to exactly 

1 m n 

v ' i=i j=i 

pairings in V(s,t). Therefore, if A is a set of multiplicity vectors, Va — {P £ T y ( s jt) I 
a(P) G A}, and A^a = {Q G A4(s,t) | a(Q) G A}, then 

lMAl ~nz^ (31) 

This holds in particular if A is the set of all nonnegative integer sequences, in which case 
V A = V(s,t) and M A = M(s,t). 

We first prove Theorem 11.31 in the case that (S 2 ,T 2 ) is insubstantial. 

Lemma 3.1. If 1 < st = o(S 2 ^ 3 ) and (S 2 ,T 2 ) is insubstantial then Theorem ] 1.3\ holds. 

Proof. Similarly to [Til Lemma 2.2], define a doublet to be to be an unordered set of 2 
parallel pairs. A double pair provides one doublet, while a triple pair provides 3 doublets. 
For the uniform probability space over V(s, t), let b r be the expectation of the number of 
sets of r doublets, for r > 0. In [TTJ Lemma 2.2] it is shown that 

bo = l, 
_ S 2 T 2 

s 3 r 3 (sf-4S 3 -2ff 2 )(r|-4r 3 -2r 2 ) 

2 2[S] 3 8[S}, 

6[6J 3 
64 = 0{sH 3 /S 2 ). 

Let p/c denote the probability that a randomly chosen pairing contains exactly k dou- 
blets, for k > 0. Then 



£M>" +1 (I) * 

r>k v 7 

and the partial sums of this series alternate above and below pk (see for example [U 



Pk 

r>k 



13 



Theorem 1.10]). Applying this, we find that 

n _ . S 2 T 2 S 3 T 3 (S 2 2 -4S 3 -2S 2 )(T 2 2 -4T 3 -2T 2 ) n( ,. 3/q2 , 

_ S 2 T 2 S 3 T 3 (S 2 -AS 3 -2S 2 ){T 2 -AT 3 -2T 2 ) 3 3 2 

Pl -2[^-2[^ qsu +o{st/s ), 

p2 = (St-,S 3 -2S 2 )(Ti-,T 3 -2T 2) + 0(s 3 t3/ , 2)) 

P3 = || + o( sW 

(The expression for p was also derived in [Tij Lemma 2.2].) The configurations defining 
these cases are, respectively, no parallel pairs, one double pair, two double pairs, and one 
triple pair. 

Applying Lemma 12.61 and (13. ip , 

M(s, t) = (1 + 0(s 3 t 3 /S 2 )) = f' ( Po + 2 Pl + Ap 2 + 6p 3 ) 

1L=i s «- llj=i t r 

(Po + 2 Pl + 4p 2 + 6p 3 + 0(s 3 t 3 /S 2 )) 



5! 



' ^ S 2 T 2 S 3 T 3 (g 2 2 -4g 3 -2g 2 )(T 2 -4T 3 -2T 2 ) A 
/ + 2[^ + 3[^ + 8[SU +0(st/S )j, 

where we have used the fact that po + 2pi + Ap 2 + 6p 3 = 1 + o(l) in the insubstantial case 
to get the second line. 

This expression is equal to the expression in Theorem 11.31 under our present assump- 
tions. (Note that since (S 2 ,T 2 ) is insubstantial, the term S , |T 2 2 /2S' 5 which appears in the 
statement of Theorem 11.31 is absorbed into the error term.) □ 

For nonnegative integers d, h, define C^ft — Cd,h( s ,t) to be the set of all pairings in 
V(s, t) with exactly d double pairs and h triple pairs, but no parallel classes of multiplicity 
greater than 3. Also define 

w(C dih )= w(P) = 2 d 6 h \C dA \. 

Pec d , h 

A special case of (13. ip . used in [TT], is that the number of 0-1 matrices in V(s, t) is 

AT"/ ,\ 1^0,0 | 

We will proceed by writing M(s, t) in terms of N(s, t), as follows. 
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Lemma 3.2. If (S 2 ,T 2 ) is substantial then 

N 2 N 3 , n 

w{C, 

d=0 h=0 



M(s, t) = N(s, t) £ £ (1 + 0(,¥/5 2 )) . 

— — W I Un n ) 

Proof. By Lemma 12.51 and ( 13. II) , 



AT 2 iV 3 



m( S , t) = 1 g 5>(<^o (i + o( S ¥/s 2 )) 

1L=i »■ llj=i j ■ d=0 ft=0 



We will evaluate the sum in Lemma 13.21 using two summation lemmas proved in [TT] 
and restated below. 

Lemma 3.3 ([TH Corollary 4.3]). Let < A\ < A 2 and B\ < B 2 be real numbers. 
Suppose that there exist integers N, K with N > 2 and < K < N , and a real number 
c> 2e such thatO < Ac < N-K + l and \BN\ < 1 for all A G [A U A 2 ] and B G [B 1 , B 2 \. 
Further suppose that there are real numbers Si, for 1 < % < N, and 7, > 0, for < % < K , 
such that Yl)=i\$j\ < Hj=QlAi\i < l f or 1 <i < N - 

Given A(i), . . . , A(N) G [Ax, A 2 ] and B(l), . . . , B(N) G [B 1 , B 2 ], define n Q , ni,...,n N by 
hq = 1 and 



rij 



for 1 < i < N , with the following interpretation: if A(i) = then rij = for i < j < N . 
Then 



N 



where 



i=0 



K 



Ih = exp(Ai ~ l 2 A 2 1 B 2 -Aj2^ A iY) - H 2e / C ) j 

3=0 

K 

E 2 = exp(A 2 - \A\B X + \A\B\ + A^ lj (?>A 2 )^ + \{2e/c) N . □ 



j=0 

Lemma 3.4 ([HI Corollary 4.5]). Let N > 2 be an integer and, for 1 < i < N, let 

real numbers A(i), B(i) be given such that A[i) > and 1 — (i — l)B(i) > 0. Define 
A 1 = mmf =1 A({), A 2 = max^^z), G x = mmf =l A{i)B{i) and C 2 = msxf =1 A{i)B{i). 
Suppose that there exists a real number c with < c < ~ such that m&x{A/N, |C|} < c 
for all A G [Ai, A 2 ], C G [Ci, C 2 ]. Define n , . . . , by n = 1 and 

JjL = M (l _ (i _ 1)BW ) 

"i— 1 b 
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for 1 < i < N, with the following interpretation: if A(i) = or 1 — (i — l)B(i) = 0, then 
Uj = for i < j < N . Then 

N 
i=0 

where 

E x =exp(A 1 - \A X C 2 ) -(2ec) N , 

E 2 = ex P (A 2 - iA 2 Ci + |A 2 C?) + (2ec) N . □ 

We obtain bounds on the ratios we require by applying results from [UJ . To begin with 
we focus on the effect of changing the number of triple pairs while keeping the number of 
double pairs fixed. 

Lemma 3.5. Suppose < d < N 2 and 1 < h < N 3 , with Cd t h 7^ 0- If {S 2 ,T 2 ) is 
substantial then 

w(C dih ) _ S 3 T 3 + 0(sH 2 (st + d + h)S) 
w(C d ,h-i) hS 3 
Proof. This follows from [HI Lemma 4.6] since, for h > 1, 

w{C dih ) _ 6 \C dth \ 
w(Cd,h-i) \Cd,h-i\ 

Next, adapting the proof of [TU Corollary 4.7] gives: 

Corollary 3.6. Suppose < d < N 2 with C dj o ^ 0- Further suppose that (S 2 ,T 2 ) is 
substantial. Then 

Proof. We will apply Lemma [3.41 Let h' be the first value of h < N 3 for which C dt h = 0, 
or h! = N 3 + 1 if there is no such value. Define a^, 1 < h < h', by 

\C dth \ _ S 3 T 3 -a h (s 2 t 2 (st + d +(h-l)S)) 



\C d ,h-i\ hS 3 

Lemma [3751 says that is bounded independently of h, d and S. 
For 1 < h < h! , define 

a/i\_ S3T3 — a h (s 2 t 2 (st + d)S) _ a h s 2 t 2 

A{ti) — , L-{n) — — — — . 

If ah < then by definition A(h) > S 3 T 3 /S 3 , and S 3 T 3 > since h < h'. Therefore 
A(h) > in this case. If ah > then C(h) > 0, which implies that A{h) > since the 
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right side of (I3.2p has the same sign as A(h) — (h — l)C(h). Therefore A(h) > whenever 
h < h'. Define B(h) = C(h)/A(h) for 1 < h < ti. Also define A(h) = B(h) = for 
h! <h< N 3 . 

Define A\, A2, C%, C 2 by taking the minimum and maximum of the A(h) and C(h) over 
1 < h < N 3 , as in Lemma I3T41 Let A 6 [Ai,y4 2 ] and C G [CijCy, and set c = j^. Since 
A = S 3 T 3 /S 3 + o(l) and C = o(l), we have that max{A/N 3 , |C|} < c for 5* sufficiently 
large, by the definition of N 3 . 

Therefore Lemma 13.41 applies and says that 

exp(^p + 0{sh\st + d)/S*)) + 0((2e/41)"») . 
h=o \ d fi\ \ / 

Finally, (26/41)^ < (2e/41) log,s < S~ 2 . Since the sum we are estimating is at least equal 
to one, this additive error term is covered by the error terms inside the exponential. This 
completes the proof. □ 

Now we must sum over pairings with no triple pairs. 

Lemma 3.7. Suppose that (S2,T 2 ) is substantial and that 1 < d < N 2 with C dj0 7^ 0- 
Then 

i^ = ^( 1 -( d - 1)B )( 1+ y 



w(C d -i,o) d 



S 2 T 2 f 1 S 2 T 2 1 : 2S 3 T 2 , 2S 2 T 3 S 3 T 3 2S 2 T 2 \ . „fsH 3 



where 

U ~^V 52 52 5 5^ ~$ 3 ~) \S 

02 J 2 J 2 '-'2 ^ 

. ^/(d-l)V (c/ - l) 2 t 2 dst(d + st) 

5 d = 0\ -2 + — 2 + — 

\ iJ 2 -L 2 o 2 l 2 

Proof. This follows from [Til Lemma 4.8] since, for d > 1, 

w(C d ,o) _ 2 \C d , \ 

Adapting the proof of [TTI Corollary 4.9] gives the following: 
Corollary 3.8. If (S 2 ,T 2 ) is substantial then 

^^w(C dh ) [S 2 T 2 S 2 T 2 „fs 3 t 3 

d=0 h=0 v u,u; v v 



□ 
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Proof. We need to apply Lemma 13.31 to the result of Lemma 13.71 and take into account 
the terms coming from the triple pairs (as given by Corollary 13.61) . 

Let d' be the first value of d < N 2 for which |Cd,o| — 0, or d' = N 2 + 1 if no such value 
of d exists. Define mo, mi, . . . , mjv 2 by 

w(Cd,o) w(C d ,h) 
md ~ w(C 0fi ) ^ w(C d , ) 

for < d < d', and rrid = for d' < d < N 2 . Then clearly 

Corollary 13.61 tells us that for d < d' we have 

m d = ^4 ex P (h5T + OisHyS 2 ) + ZdsHysA (3.3) 
w{C ,o) V s J 

where £o = and in general ^ = 0(d). (Note that (13.31) is also true for d' < d < N 2 , 
since both sides equal zero.) If a is a constant such that < ocd for < d < d', then 



^ + 0(sH*/S 2 ) £ n,(-l) <^m d <exp(^ + 0(sV/5 2 ) £ n d {l) 

' d=0 d=0 ^ ' d=0 

(3.4) 

where 

^d(a^) = — 77T~l exp(xads 2 t 2 1 S 2 ) . 
w{C 0: o) 

Next we note that, for x G {—1, 1}, n (x) = 1, and for 1 < d < d', 

" A ' r) 2A(d)(l-(d-l)B)(l + 5 d ) 



n d -i(x) 

with A(d), B, and 5d satisfying the expressions given in the statement of Lemma 13.71 
This follows since the factor exp(xas 2 t 2 / S 2 ) is covered by the error term on A(d). For 
d! <d<N 2 define A(d) = 0. 

Now let A\ = A\(x) = min.d2A(d), A 2 = A 2 (x) = m&Xd2A(d), where the maximum 
and minimum are taken over 1 < d < N 2 . Also let B\ = B 2 = B, and K = 3, and define 
c = S* 1 / 4 if S 2 T 2 < S 7 / A and c = 41 otherwise. The conditions of Lemma [3.31 now hold as 
we will show. Let A G [Ai,^] be arbitrary. 

Clearly c > 2e. If S 2 T 2 < S 7 / 4 then N 2 = 22. Using the condition S 2 T 2 > (st)^ 2 S 
implied by the substantiality of (S 2 ,T 2 ), we find that Ac = 1 + o(l). For S 2 T 2 > S 7 ^ 4 , 
Ac = 4AS 2 T 2 / S 2 {\ + o(l)). It is also easy to check that BN 2 = o(l). Thus, in all cases 
we have that Ac < N 2 — 2 and \BN 2 \ < 1 for sufficiently large S. 
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If d = 0(S 2 T 2 /S 2 ) then 

t 2 &T„ «rf<7?7? c 2,2c^ T „\ / c 3,3\ 

= 0(1), 



S 2 



while if d < [log S 1 ] then 

^| = ^^_ + _^ + ^_ + ^_) =0(1 ). 



=1 



Finally, for 1 < < A^, we have 



EN^(^-i)^ + ;)) + o(t||) + o(E rfs2(2 



S2T2 

/- ^ />(/>'+ l).s/ £;(£; + l)s 2 t 2 



5*2^2 S2T2 



K 



J'=0 

where 



s 2 t 2 \ „ / s 2 t 2 s 2 t 2 \ / ' s 2 t 2 st 



7o = 0, 7l = — - , 72 = — + — + — - , 7 3 = 0-.+ 



C T 1 / ' ' \ C2 T->2 C / ' ' \ C2 T^2 C T" 1 

02 I 2/ \>->2 l 2 ^>2 1 2 / \^>2 l 2 ^>2 1 2, 

Since A^|(s 2 / '5f + t 2 /T 2 2 + st/S2T 2 ) = o(l), which is easily checked, it follows that 
J^lo 7j [&] j < 1/5 for 1 < k < N 2 , when S 1 is large enough. 

Therefore the conditions of Lemma [3731 hold, and we conclude that each of the bounds 
given by that lemma for Y2d=o n d( x ) has the form 

exp (a - \ A 2 B + O [a 3 B 2 + £ -yj(3A)A )+€> ({2e/c) N *), 

where A is either Ai or A 2 . A somewhat tedious check shows that 

3 

0(A 3 B 2 ) +J2lj^ A ) j = 0(sH 3 /S 2 ). 
3=0 

Next consider the error term 0((2e/c) 7V2 ). If N 2 = 22 then (2e/c) N2 = ^eS^ 1 / 4 ) 22 = 
0(S- 2 ), while in the other cases we have (2e/c) N2 = (2e/Al) N2 < (2e/41) log5 = 0(S~ 2 ). 
Since no = 1, this additive error term is covered by a relative error of the same form. 
Therefore, each of the bounds on Y2d=o n d( x ) nas the form 

e Xp (^I^ + (f))=e X p(fi + fi-^ +0 (^ 
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Modulo the given error terms, the final expression does not depend on x, nor on whether 
we are taking a lower bound or upper bound in Lemma [3. 31 To complete the proof, just 
apply (1331). □ 



Corollary 13.81 and Lemma I3~2l together prove Theorem 11.31 in the substantial case. The 
insubstantial case was already proved in Lemma 13. 11 



4 Alternative formulation 

We now derive an alternative formulation of Theorem 11.31 Recall the definition of and 
z>fc given in the Introduction. 

Corollary 4.1. Under the conditions of Theorem \1.3[ 

n+Si— 1\ -pT f m+tj — 1 



n " T r n 



Si / - LJ r V t 



M (s,t) = — 



mn+S—1 
S 



x exp (1 — /x 2 )(l - 2 ) - 



1 3 - fi 2 i> 2 \ 



2 AS J 
[1 - /2 2 )(3 + £2 - 2fi 2 u 2 ) (1 - i>2)(3 + u 2 - 2p, 2 i>2) 



An Am 
+ (l-3/j 2 + 2^ 3 )(l-3z> 2 + 2i> 3 ) + Q fsH 3 



12S V S 

Proof. By Stirling's formula or otherwise, 

iV+a;— 1\ N* f[x] 2 [x] 3 [x} 2 4/ 



x J d eXP l^-o^-4^ + ^ 
as N — > 00, provided that the error term is bounded. This gives us the approximations 
«/ n+a ._l\ _n^_ [S 2 _S 2 _ A.n^ 

" / m +tj-l\ _rrf_ (T^ Tj T^ + nf^f. 

V *i / n^i !eXP V2m 4m 2 6m 2 V ^ 2 

mn+S-l\ (mn) s ( S 2 S S 3 ^fsH 3 

ex P n n 7^T^ + ° 



S J S\ \2mn 2mn Qm 2 n 2 \ S : 

Substitute these expressions into Theorem 11.31 and replace S 2 , S 3 ,T 2 ,T 3 by their equiva- 
lents in terms of jj 2 , fa, v 2) £3. The desired result is obtained. □ 
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As noted in the Introduction, Theorem 11.31 establishes the conjecture recalled after 
Theorem 11.11 in some cases. Using Corollary 14.11 the following is easily seen. (Note that 
A*2 — A*3 — ^2 = ^3 = in the semiregular case.) 

Corollary 4.2. If s = s(m,n) and t = t(m,n) satisfy ms = nt and st = o((mn) 1 '' 5 ), 
then 

A(m,s;n,t) = ^±% + (l)). □ 

Most of the terms inside the exponential of Corollary 14.11 are tiny unless at least one 
of fa, &2 is quite large (that is, the graph is very far from semiregular). In particular we 
can now prove Corollary 11.51 which was stated in the Introduction. 

Proof of Corollary \1.5l It is only necessary to check that the additional terms in Corol- 
lary HJ] have the required size. It helps to realise that /2 2 < s, [fa] < sfM 2 , 2 < t and 
\h\ < tD-2. □ 



A random nonnegative mxn matrix with entries summing to S is just a random com- 
position of S into mn parts. (A composition is an ordered sum of nonnegative numbers.) 
In particular, for 1 < i < m the row sum Sj satisfies 

Pr(s , = »)=(* + I - i )( S - k + 1 >" - !) /( S + ™" - l ) ( <k<S). 

From this we can compute the following expected values. 

^ n(m — 1) ^ A m(n — 1) 



mn +1 " mn + 1 

n(m — l)(m — 2){mn + 25*) m(n — l)(n — 2)(mn + 25") 

E ^ 3 = ? r^y\ ' Ez/ 3 = 7 — — 77 • 

m[mn + 1) (mn + 2) n[mn + lj(mn + 2) 

The first two expectations suggest that the argument of the exponential in Corollary 11.51 
is close to with high probability for such a random matrix. We will prove this in a future 
paper, and note that the result gives a model for the row and column sums of random 
matrices. 



5 Restricted sets of allowed entries 

Given a subset J of the nonnegative integers, let M(s, t, J) denote the set of matrices in 
M(s, t) with all entries in the set J. Let M(s, i, J) = \M(s, t,J)\. By generalising the 
techniques of the preceding sections, we can find an asymptotic expression for M(s, t, J) 
whenever 0, 1 G J . 
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Lemma 5.1. Let J C N with 0,1 e J. Define X2 = if 2 J, \2 = 1 »/ 2 G J", and 

similarly X3- Then 



S 2 



S 3 



S 3 



3+3 



M(s, t, J) = AT( S) t) exp ( + X2 S 4± + ( X3 - X 2 )^5 + O / 



S 2 



5! 



GXP (X2-2)-55- + (X2-2) _ OT" + (X3 - X2 + 



2 S 2 

"2) 



S 3 



S 2 T 2 (5*2 + T 2 ) S 2 T 3 -\- S 3 T 2 S 2 T 2 



AS 4 



2S A 



+ 



2S 5 



37 S 3 
s 3 ^ 



S 2 



Proof. Our general approach will be similar to that we used for Theorem 11.31 but the 
methods of Section [2] will need significant modification. The source of the problem is that 
a D-switching may introduce an entry that is not in J . 

For Q E M.(s, t) and i > 1, let ni(Q) be the number of entries of Q equal to i. Also 
let n> 5 (Q) = J2 i>5 ni(Q). Define N 2 and N 3 as in Section [2] when (S 2 ,T 2 ) is substantial, 
and N 2 = 2 and N 3 = 1 when (S 2 , T 2 ) is insubstantial. Consider the following subsets 
of M(s,t): 

M + = M(s,t, t 7U{4,5,6,...}), 

.M = M(s,t,J), 
M- = {QeM + \ n 2 (Q), n 3 (Q), n 4 (Q), n> 5 (Q) < S 5 / 6 }, 
M* = {Qe M(s,t, J n {0,1,2,3}) I n 2 (Q) < N 2 , n 3 (Q)<N 3 }. 

Also define the cardinalities M + , M, M~ , M* , respectively. By monotonicity, we have 
M* < M < M ■ and M* < M~ < M+. 

We now employ switchings to establish that M + -M~ = o(M~—M*) and M~-M* = 
0(sH 3 /S 2 )M*, from which it follows that M = (l + 0(sH 3 /S 2 ))M*. 

We start with the claim that M + —M~ = o(M~ —M*). This is accomplished with a 
succession of four types of switchings, illustrated by the following operations on subma- 
trices. 



qi 

q 2 

A o\ 

4 

4 

\0 4/ 



1 

1 92-1 

(1 1 1 i\ 
1111 
1111 
1 1 1 i) 



Qi, <?2 > 5 
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3 0\ fl 1 1\ 

3 i — ^ I 1 11] 
3/ \1 11/ 




To be precise, given any Q £ 7W + \.M~, we apply the first switching until n> 5 (Q) < S 5 / 6 , 
then the second switching until n^Q) < S 5 ^ 6 , then the third switching until ris(Q) < S 5 ^ 6 , 
then finally the fourth switching until n2(Q) < S 5 ^ 6 . Since none of the switchings can 
undo the work of a previous switching, the end result is a matrix in M~ \A4*. (Note 
that in the resulting matrix R, at least one of n 2 (R), n 3 (R), ra 4 (i?), n> 5 (R) has value 
55/6 _ 0(1). This implies that R g" M* .) 

Quantitatively, in each of the four phases, we remove a part of Ai + \ M.~ which is 
small compared to Ai~ \ Ai*. As a representative example, take the second switching. 
By counting similarly to Lemma [2T2| the second switching can be applied to Q in at least 
(n 4 (Q) — O(st)) 4 ways, and the inverse can be applied in at most Ss 3 t 3 ways. For n 4 (Q) > 
S" 5 / 6 , the condition s 3 t 3 = o(S 2 ) implies that Ss 3 t 3 = o[(n 4 (Q) — O(st)) 4 ), so the ratios 
that play the part denoted by a(t> j) in Theorem 12. II are all o(l). The same type of 
argument applies in all four phases, which establishes the bound M + —M~ = o(M~ —M*). 

When Q 6 M.~\J\A* , we know that rii(Q) = S — o(S). We can now continue precisely 
as in Lemmas 12. 3\\2. 51 using D-switchings restricted to gi = • • • = — 1. This restriction 
ensures that D-switchings only create entries with value equal to or 1. The various 
switching counts can be taken as essentially the same as before, since all but a vanishing 
fraction of the non-zero entries are 1. We conclude that M~ —M* = 0(s 3 t 3 / S 2 )M* which, 
as noted above, implies that M = (l + 0(s 3 t 3 /S 2 ))M*. 

Having now reduced the task to evaluation of M*, we can complete the proof following 
Lemma [3TT1 in the insubstantial case, and Section 4 in the substantial case. In Lemma [3TT1 
the only modification is to replace the expression p + 2pi + 4p 2 + Qp3 by po + 2\2 Pi + 
4X2 P2 + 6x3 P3- 

Now suppose that (5*2, T 2 ) is substantial. If X2 = X3 then the result is given by either 
Theorem 11.21 or Theorem 11.31 If X2 = and X3 = 1 then the result follows from applying 
Corollary 13.61 with d = 0, since arguing as in Lemma [3.21 shows that 

M(s, t, J) = N(s, t) £ ^((1 + 0(s 3 t 3 /S 2 )) 
in this case. Finally, if X2 — 1 and X3 — then 

M(s, t, J) = N(s, t) ^4 (1 + 0(sh 3 /S 2 )) 
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so in place of (13.31) we simply have m d = u>(C d)0 )/u>(Co,o) — n d(fy for < <i < N 2 . 
The remainder of the proof is identical except that there is no need to apply (13. 4p at 
the end. □ 
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